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ON THE PERIODIC SOLUTIONS OF A 
GENERALIZED SMOOTH OR NON-SMOOTH 
PERTURBED PLANAR DOUBLE PENDULUM 

JAUME LLIBRE 1 , DOUGLAS D. NOVAES 2 AND MARCO ANTONIO TEIXEIRA 2 



Abstract. We provide sufficient conditions for the existence of periodic solutions with small 
amplitude of the non-linear planar double pendulum perturbed by smooth or non-smooth 
functions. 



£/} \ 1. Introduction and statement of the main results 



q 

C| ' distance between a point P (called pivot) and mi and the distance between mi and m 2 are 



We consider a system of two point masses mi and 7712 moving in a fixed plane, in which the 



fixed, and equal to l\ and l 2 respectively. We assume the masses do not interact. We allow 
gravity to act on the masses mi and m 2 . This system is called the planar double pendulum. 

The position of the double pendulum is determined by the two angles 0i and 02 shown in 
Figure [T] The corresponding Lagrange equations of motion are 

CN '. 

k*" ■ (mi + m 2 )h4>i + m 2 l 2 4> 2 cos(0i - 2 ) + (mi + m 2 )gsm(4>i) 

00 ■ 
0\ ; (!) +m 2 l 2 <j>l sin 0i - 4>2) = 0, 

Q \ m 2 h4>i cos(0i - 4> 2 ) + m 2 l 2 4> 2 + m 2 gsm((j)2) + m 2 li4>\ sin(0i - (j> 2 ) = 0, 

where g is the acceleration of the gravity. For more details on these equations of motion see [S] . 
Here the dot denotes derivative with respect to the time T. 

The authors in [8] have studied, in the vicinity of the equilibrium 4>i ~ 4> 2 — 0, the persistence 
of periodic solutions of system (TTJ) perturbed smoothly in the particular case when mi = m 2 
and /1 = l 2 . Now mi,m 2l li and l 2 can take arbitrary positive values and we shall study the 
S^ \ periodic orbits of system (TTJ) which persist under smooth and non-smooth perturbations. 

Denote the expressions for <fii and 4> 2 in (TTJ) respectively by Hi (0i ,<j>i,</>2,<j>2) an d H 2 (</>i ,<j>i,4>2 
In this paper we shall consider the perturbed problem 



(2) 



ffi(0i,0i,</> 2 ,</> 2 ) +£ ( A (t,(f>i ,4>i,<f>2, <fo) + Fa(t, 0i, 0i,02, 2 )sgn((/)i)j 
+e 2 (i?i(t, 0i, 0i, 02, 2 ) + R 2 {t, 0i, 0\, 2 , 2 )sgn(0i)) + C(e 3 ), 

•^2(01,01,02,02) +e(P 3 (t, 01, 01, 02,02) + F±(t, 01, 01,02, 2 )sgn(0 

+e 2 ^R 3 (t, 0i, 0i, 02, 2 ) + & 4 (*, 0i, 0i, 02, 2 )sgn(0 2 )) + 0{e 3 ). 
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FIGURE 1. The planar double pendulum. 

The function sgn(z) denotes the sign function, i.e. 

-1 ifz<0, 

sgn(z) = { if z = 0, 

1 if z > 0. 

Here the smooth functions Fi and Ri for i = 1, 2, 3, 4 define the perturbation. These functions 
are respectively Tp i -periodic and T Ri -periodic in t and respectively in resonance PFi'-QFi and 
PRi'-QRi with some of the periodic solutions of the linearized unperturbed double pendulum, 
being p and q relatively prime positive integers for p = pp i , p Ri , q = qp t , q Ri and i = 1,2,3,4. 
We also assume that F t (t, 0, 0, 0, 0) = for i = 1,2, 3, 4. 

Remark 1. For simplicity, we can assume that the functions Fi and Ri for i = 1,2,3,4 
are T -periodic with T = pTj for some integer p where Tj for j = 1,2 are the periods of the 
solutions of the linearized unperturbed double pendulum. Indeed, if we take p the least common 
multiple among pp i and p^ for i = 1,2, 3, 4, then there exists integers np i and UR t such that 
P = np i pp i = np ii pp ii for i = 1, 2, 3, 4. Hence 

rr, PF X „ PR, rj, 

pTj = n Fi q Fz T,j = n Ri q Ri Tj. 

QF, qR z 

Fori = 1,2,3,4, and j = 1,2. 

Note that the functions Fi and Ri for i = 1,2, 3, 4, can be taken in a certain way arbitrary, 
i.e., only assuming some hypotheses. It makes us able to provide, in a physical context, real 
meaning for these functions. In our case, since we are working with discontinuity in the variables 
6\ and 02, the functions F\, F^, R\ and R2 could model the escapement for the particle mi, and 
the functions -F3, F4, R3 and R4 could model the escapement for the particle 7712. If we work 
with discontinuity in the variables 6[ and 9' 2 , instead with discontinuity in the variables B\ and 
6*2, the respective functions could model the Coulomb Friction. We also can work composing 
these two phenomena. For more details on physical systems with discontinuous models see, for 
instance, Q] and [2]. 

Now, we follow the steps: 



(i) proceed with the change of variable <pi = e9\ and (f> 2 = s9 2 ; 

6i = -HrfeOi, eOi, s0 2 , s0 2 ) + A (t, eOi, edi, £02, efe) + F 2 (t, eO 1 ,e8 1 ,e0 2 , s6 2 )sgn(6i) 

£ 

+e (&! [t, s9 1 ,e8 1 ,s9 2 , e9 2 ) + R 2 (t, e6 x ,e6\,e6 2 , e0 2 )sgn(0i)) + 0(e 3 ), 
k= -H 2 {e9 ll e9 1 ,e9 2 ,e9 2 ) + F 3 {t,s9 ll s9 1 ,e9 2 ,e9 2 ) + F 4 (t,Edi,ed u e9 2 ,ed 2 )sgn(6 2 ) 

_3\ 



+e [Ri (t, e9 1 ,e9 1 ,e9 2 , ed 2 ) + R 2 (t, e$i, s6 x , e6 2 , e0 2 )sgn(9 2 )\ + 0(e 

(ii) expand in Taylor series, for e = 0, the expressions of 0\ and 2 ; 
(iii) Take a new time t given by the rescaling t = ar, with a = ^/l\ mi/(gm 2 ); 
(iv) and finally, denote 

mi+m 2 h(mi+m 2 ) 

a = > 1 and b = > 0. 

m 2 l 2 m 2 

Hence, we obtain the following equations of motion for the double pendulum 

(3) 

0" = -O0i + 9 2 +e(K 1 (r) + F 1 (t,9 1 ,6' 1 ,0 2 ,9' 2 ) + (A' 2 (r) + F 2 (T,9 1 ,9[,9 2 ,9 2 ))sgn(9 1 )) 

+e 2 R 1 {T,9 1 ,9' l ,9 2 ,9' 2 ,e), 

9'i = 60! -b9 2 + e (K 3 (t) + F 3 (t, 0i,9[,6 2 ,9' 2 ) + (K 4 (r) + F 4 (t, 6i,9[,d 2 , 9' 2 )) sgn(fe)) 
+£ R 2 (t,9i,9 1 ,9 2 ,9 2 ,£), 

where now the prime denotes derivative with respect to the new time t. Here the functions Fi, 
for i = 1, 2, 3, 4, are linear in the spatial variables, and are given by 

F i (T,9 1 ,0' 1 ,0 2 ,9' 2 ) = d\{T)9 1 + dUT)9' l +d l {Tf9 2 + d l {T)%. 

with 

f)F ^ 8F 

4(r) = a 2 -^(ar,0), d 2 (r) = a-f (or.O), 

dUr)=a 2 -^(ar 7 0), 4(r) = a-f (ar,0), 

<7</>2 <9</> 2 

and A'((t) = a 2 A l (ar, 0, 0, 0). Observe that, for i,j = 1,2,3,4, d|(r) and Aj(r) arc T/a- 
periodic functions. 

The objective of this paper is to provide a system of equations whose simple zeros provide 
periodic solutions (see Figure [2]) of the perturbed planar double pendulum ([2]). 

In order to present our results wc need some preliminary definitions and notations. 

The unperturbed system Q has a unique singular point, the origin with eigenvalues ±wi i, ±cu 2 i, 
where 

\/a + b - a/A \fa + b + \/A 

Ul = 7=2 ' "* = 7=2 ' 

with A = (a — b) 2 +Ab > 0. Consequently this system in the phase space (#i , 6[ , 6 2 , 9' 2 ) has two 
planes filled with periodic solutions except the origin. The periods of such periodic orbits are 

2tt 2tt 

T\ = — or T 2 = —. 

Ul\ L0 2 

These periodic orbits live into the planes associated to the eigenvectors with eigenvalues ±cji i 
or ±uj 2 i, respectively We shall study which of these periodic solutions persist for the perturbed 
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FIGURE 2. Periodic solution of the perturbed system ([2| converging to the origin, 
when e — > 0. 



system (|2|) when the parameter e is sufficiently small and the functions of perturbation Fi and 
Ri for i = 1, 2, 3,4 have period either paT\, or paT2, with p positive integer. 

Remark 2. We say that the Crossing Hypothesis is satisfied if there exists a compact set 
Del such that every orbit starting in D reaches the set of discontinuity only at its crossing 
regions (see Appendix A). 

Let Xx ,Yo( T ) be the periodic function 

Xx ,Y (t) = Y cos (wi r) + X sin (wi r) , 

then we define the non-smooth function J-"i(Xo,lo) by 

rpTi 

in ( Wl r) (2&(F X + A^(r)) + (F 3 + K 3 (t)) (a - b + \/Aj J dr 

- / sin(wiT)(26( J F 2 +A 2 (T)) + (F4 + A4(r))( a -fe + VA)) sgn(Xx ^ (r)) 

and the non-smooth function J-2(Xo,lo) by 
rpTi 



(4) " rP T; 



dr. 



cos (wi r) (^26(Fi + K x (t)) + (F 3 + K 3 (t)) i y a-b + VDjj dr 

(5) J ° r pTl / / xx 

+ / cos (lj! t) (2b(F 2 + K 2 (t)) + (F 4 + K 4 (t)) (a - b + Vd)) sgn(X Xo , y (t)) dr. 



where 

(0) 



F i = F i (j,Ai,B 1 ,C 1 ,D 1 ) 



for i = 1,2,3,4 with 



A 1 = - — (X cos (wi r) + y sin (cji r)) , 

Z LO\ 

(-a + b + y/A\ 

Bi = — (Y cos (uji t) - X sin (w x r)) , 

Ci = — (Xo cos (wi r) + y sin (wi r)) , 

Z?i = Fo cos (wi r) — X sin (wi r) . 
A zero (Xq, F *) of the system of the non-smooth functions 
(7) F 1 (X ,Y ) = 0, T 2 (X ,Y )=0, 

such that 



dct 



(x ,y )=(x;,y*)J 



K d(X ,Y ) 

is called a simple zero of system J7} . 

Our main result on the periodic solutions of the non-smooth perturbed double pendulum 
@ which bifurcate from the periodic solutions of the unperturbed double pendulum ((T|) with 
period T\ traveled p times is the following. 

Theorem 1. Assume that the functions Fi and Ri of the non- smooth perturbed double pendulum 
(|2|) are periodic in t of period paT± with p positive integer. Also assume that the Crossing 
Hypothesis (see Remarket is satisfied. Then for \e\ > sufficiently small and for every simple 
zero (Xq,Yq) ^ (0,0) of the non-smooth system JjJ) such that the orbits pass by D, the non- 
smooth perturbed double pendulum |(5J) has a paT\ -periodic solution ((f>i(t,e),4>2(t,£)) — > (0,0) 
when e — > 0. 

Theorem [T] is proved in section [2J Its proof is based in the averaging theory for computing 
periodic solutions, see the Appendix B. 

Note that the periodic solution given in Theorem [T] is a periodic solution bifurcating at e = 
from the equilibrium of system @ localized at the origin of coordinates. For |e| > sufficiently 
small this orbits is close to the plane defined by the eigenvectors of the eigenvalues ±zwi. 

We provide an application of Theorem [T] in the following corollary, which will be proved in 
section [3) 

Corollary 2. Suppose that Fi = y/a 2 +fi, F% = w/a 4-/3, and fi, F2, fz, and F4 has no linear 
term. Also suppose that R2 = l/a 2 + r2, R4 = l/a 2 + 7"4 and R\, r-i, R3, and r^ has no constant 
term. Moreover, assume that all functions are paTi/q -periodic. Then the differential system 
(J2|) for \e\ > sufficiently small has one paTi -periodic solution {4>i(t,e),4>2{t,e)) — > (0,0) when 
e^0. 

Now let Z Zo ' Wo (t) be the periodic function 

z z ,Wo ( T ) = Wo cos (^ r ) + z a sin (wa t) , 
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then we define the non-smooth function J rl (Zo, Wo) by 
(8) 



f 2 sin(w 2 T)(-2fe(F 1 +A" 1 (T)) + ( J F H 3 + A' 3 (T))(-a + 6 + VA)) dr 

sin (wj, r) (2b(F 2 + K 2 {t)) + (F 4 + A' 4 (r)) (-a + 6 + \/a)) sgn(Z z °- w ° (r)) dr, 



pT 2 





and the non-smooth function J r2 (Zo, Wo) by 

pT-2 , / \\ 

cos (W2 r) f-26(F! + A^r)) + (F 3 + K 3 (t)) (a - b + y/DJJ dr 

,pT 2 , _ _ , xx 

+ / cos ( Wa r) [2b(F 2 + K 2 {t)) + (A 4 + A 4 (r)) [a-b+ y/Djj sgn(Z z °> w ° (r)) dr. 



(9) •"> 



'o 
where 

for i = 1,2,3,4 with 



F i = F i (t,A 2 ,B 2 ,C 2 ,D 2 ) 



a-b+ V A 
A 2 = -+ — '- (Z cos (w 2 i) + Wo sin (w 2 £)) , 

Z 6^2 



-6 + a/A 



a 
B 2 = — ^ — '- (Wo cos (W2 *) - ^o sin (w 2 *)) , 

C 2 = — (Zq cos (w 2 t) + Wo sin (u 2 t)) , 

oj 2 

D 2 = Wo cos (lu 2 t) — Zq sin (w 2 £) ■ 

Consider the non-linear and non-smooth system 

(10) F 1 {Z ,Wo)=0, -F 2 (Z ,Wo)=0. 

Our main result on the periodic solutions of the non-smooth perturbed double pendulum 
@ which bifurcate from the periodic solutions of the unperturbed double pendulum (|3|) with 
period T 2 traveled p times is the following. 

Theorem 3. Assume that the functions Fi and Ri of the non- smooth perturbed double pendulum 
@ are periodic in t of period paT 2 with p positive integer. Also assume that the Crossing 
Hypothesis (see Remarket is satisfied. Then for e > sufficiently small and for every simple 
zero (Zq, Wq) ^ (0,0) of the non-smooth system (|10p such that the orbits pass by D, the non- 
smooth perturbed double pendulum IJ2J has a paT 2 -periodic solution (<f>i(t,e),<f) 2 (t,e)) — > (0,0) 
when e — > 0. 

Theorem [3] is also proved in section [2] 

Again the periodic solution given in Theorem [3] is a periodic solution bifurcating at e = 
from the equilibrium of system $Z§ localized at the origin of coordinates. For \e\ > sufficiently 
small this orbits is close to the plane defined by the eigenvectors of the eigenvalues ±zu; 2 . 

We provide an application of Theorem [3] in the following corollary, which will be proved in 
section [3l 

Corollary 4. Suppose that F\ = y/a 2 +fi, A3 = w/a 2 +fs, and fi, F 2 , fa, and A 4 has no linear 
term. Also suppose that R 2 = l/a 2 + r 2; A 4 = l/a 2 + r 4 and Ri, r 2 , R3, and r 4 has no constant 
term. Moreover, assume that all functions are paT 2 /q -periodic. Then the differential system 



(J2J) for \e\ > sufficiently small has one paT 2 -periodic solution (</>i(£,e), <fi 2 (t,£)) — > (0,0) when 

2. Proofs of Theorems [1] and [3] 

Introducing the variables (x,y,z,w) = (9i,9' 1 ,9 2 ,0' 2 ) we wr ite the differential system of the 
non-smooth perturbed double pendulum (|3]) as a first -order differential system defined in IR 4 . 
Thus we have the differential system 



x = y, 

1/ = -ax + z + e (Ki(t) +F 1 (T,x,y,z,w) + (K 2 (t) + F 2 (r,x,y, z,w))sgn(x)) 

+e 2 Ri(r,x,y,z,w,e), 
z' = W, 
w' = bx- bz + e (K 3 (t) + F 3 (T,x,y,z,w) + (K±(t) + Fi{r,x,y, z, w))sgn(z)) 

+£ 2 R 2 (t, x, y, z, w,£). 



(11) 



System (jTTJ) with e = is equivalent to the unperturbed double pendulum system ([3]) , called 
in what follows simply by the unperturbed system. Otherwise we have the perturbed system. 

Instead of working with the discontinuous differential system (|IIj) we shall work with the 
smooth differential system 

x' = y, 

1/ = -ax + z + £(Ki(t) +Fi(r,x,y,z,w) + (K 2 (t) + F 2 (r,x,y, z,w)) s s (y)) 

+£R 1 (T,x,y,z,w,£), 
(LZ) 

z = w, 

w' = bx- bz + £ (K 3 (t) + F 3 (r,x,y,z,w) + (A" 4 (t) + F 4 (r, x,y, z, w)) s$(w)) 

+£R 2 (r,x,y,z,w,£). 

where s$(x) is the smooth function defined in Figure [51 such that 

lim s$(x) = sgn(x). 
6— >0 



sign(x) 



1 


. 


s s (x) 


-S 


f 




6 


x 



FIGURE 3. The functions sign(a;) and ss(x). 
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We shall write system ([12]) in such a way that the linear part at the origin of the unper- 
turbed system will be in its real normal Jordan form. Then, doing the change of variables 
(r, x, y, z, w) — > (r, X, Y, Z, W) given by 



/ 



(13) 



1 x" 

Y 
Z 



\ 



bui\ 




b 



b 


Co 

UJ'2 


1 ffl-Ot V^A 1 


n 


Va 



bu> 2 


2VA 


(-a + b + Va) 


a-6+ VA 
2VA 

n 







2VA 



-a + 6 + VA 



\ 



( , \ 



X 

y 

z 



V w J 



2\/A / 



the differential system (|12l) becomes 

x' = Wl y, 



(14) 



r 



z' 



1 



1 



'2VA 



26 ^(r) + Fi + {K 2 (r) + F 2 ) s s (A) 



2VA 

U) 2 W, 



a-b + VA) [K 3 (t) + F 3 + (K 4 {t) + Fi) s s (C)) +e 2 i?i 



W' = -w 2 Z + e 



2VA 



-2MAi(t) + F 1 + (K 2 (t) + F 2 ) s s (A) 



+£ ^- a + b 



K 3 (t) + F 3 + (K 4 (t) + Fi) ss(C) + e'Ra, 



where Fi(t,X,Y,Z,W) = Fi(t,A,B,C,V) for i = 1,2,3,4 with 



A 



-a+b+ VA 
2buj 1 



X 



-& + VA 
2buj 2 



-a + 6+VA) (a-6+VA) 
6 = ^ '-Y - -^ '-w, 



1 



26 

1 



26 



-X + —Z, 

Wl W2 



p = y + w, 

and i?i = Ri(X,Y, Z,W,e). Note that the linear part of the differential system (TT4")) at the 
origin is in its real normal Jordan form. 

Lemma 5. The periodic solutions of the differential system (|14[) with e = are 

■3txo,Yb (r) = X cos (cVi r) + Y a sin (wi r) , 

yjfo.yb ( r ) = *o cos (wi r) - X sin (wi t) , 
(15) 

^x ,r ( T ) = °> 

Wx o ,y o (t)=0, 



of period T\ , and 



Y z °' w °(t) = 0, 
(16) V ' 

Z Z ,W ^ = Zq cog (- W2 T ) + ^ gin (- W2 T ) ) 
^Zo.Wo ( T ) = ^ Q cos ( W2 T ) _ ^ Q sin ( W2 T ) ; 

0/ period T 2 ■ 

Proof. Since system (jT^J) with e = is a linear differential system, the proof follows easily. □ 

Proof of Theorem [7J Assume that the functions Fi and Ri of the non-smooth perturbed double 
pendulum with equations of motion ((3]) are periodic in t of period paT\ with p positive integer. 
Thus Ki and Fi are pT\ -periodic functions, i.e., the differential system ([3]) and the periodic 
solutions (|T5|) have the same period pT\. 

It is well known that a Poincare map defined in a smooth differential system is smooth. So 
the Poincare maps associated to the periodic orbits of the differential system (|T2"j) arc smooth. 

The Poincare maps, restricted at D, associated to the periodic solutions of the non-smooth 
differential system (jTTJ) , which are perturbations of the periodic solutions (fT5|) are also smooth. 
Since the orbits starting in D reaches the discontinuity set only at the of crossing region (see 
Appendix A), such Poincare maps are compositions of smooth Poincare maps. In a similar way 
it follows that the Poincare maps, restricted at D, associated to the periodic solutions of the 
non-smooth differential system (fTTj) , which are perturbations of the periodic solutions (TT6"|) are 
also smooth. 

We can use Theorem [6] (see the Appendix B) for computing some periodic solutions of the 
smooth systems. The periodic solutions are zeros of the displacement function, which is the 
Poincare map associated to periodic solutions minus the identity. In fact, the non-linear function 
(|25[) whose zeros can provide periodic solutions, is the first term of order e of the displacement 
function. See for more details the proof of Theorem [6] in [3] . 

Since the Poincare maps associated to periodic solutions of system (|11|). coming from the 
perturbed periodic solutions (fT5|) or (|T6|) . are smooth and these Poincare maps are the limit 
of the Poincare maps associated to the smooth system (|12l) , for which we can use Theorem [51 
it follows that we also can use Theorem [5] for computing some of the periodic solutions of the 
non-smooth system (fTTj) . In other words, we can apply Theorem [6] to the smooth systems (|12[) 
and then pass to the limit, when 5 — s- 0, the function (|25p for obtaining a function whose zeros 
can give periodic solutions of the non-smooth system (jlip . 

We shall apply Theorem [5] of the Appendix B to differential system (fTT| . We note that 
system (fTT)) can be written as system ((22J) taking 



' A N 

Y 
Z 



I 



, t = r, G (i,x) 



uiY, 



\ 



-wiX, 

ui 2 W, 

y -oj 2 z j 
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Gl(i,x) = 



-= ( K^t) + F t + (K 2 (t) + F 2 ) ss(A) 



K 3 (t) +F 3 + (Iu(r) + Ft) s 5 (C) 




2VA 



--= (#i(r) + Fi + (K 2 (t) + F 2 ) s s {A) 



and G 2 (t,x, e) 



( ° ) 


Ri 





UJ 



We shall study which periodic solutions (fT5|) of the unperturbed system (|14[) with e = can 
be continued to periodic solutions of the perturbed system (TH| for e^O sufficiently small. 

We shall describe the different elements which appear in the statement of Theorem [6] in the 
particular case of the differential system (jT^J). Thus we have that f2 = R 4 , k — 2 and n = 4. Let 
ri > be arbitrarily small and let r 2 > be arbitrarily large. We take the open and bounded 
subset V of the plane Z = W = as 



v = {(x ,y ,o,o) g R 4 : n < y^ 2 + r 2 < r 2 }. 

As usual C1(V) denotes the closure of V. If a = (Xo, Yj), then we can identify V with the set 

{a G R 2 : n < ||a|| < r 2 }, 

here || ■ || denotes the Euclidean norm of R 2 . The function /3 : C1(V) ->■ R 2 is /3(a) = (0,0). 
Therefore, in our case the set 

Z = {z a = {a, (3(a)) , a G C1(F)} - {(*„, F , 0, 0) G I 



ri < V^o + ^ < M- 

Clearly for each z a <G Z we can consider the periodic solution x(r, z Q ) = (X(t),F(t), 0,0) 
given by l[T5]) with period pT\ . 

Computing the fundamental matrix M Za (r) of the linear differential system ([T4")| with e = 
associated to the T-periodic solution z a = (Xq,Yo,0,Q) such that M Za (0) be the identity of 
R 4 , we get that M{t) = M Za (r) is equal to 



/ 



cos (u)\ t) sin (bj\ t) 

- sin (wi r) cos (wi r) 







\ 



cos (W2 t) sin (cj2 r) 

— sin (w2 r) cos (W2 t) j 
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Note that the matrix M Za {r) does not depend on the particular periodic solution x(r, z a ). 
Since the matrix 



/ 



M- 1 (0)-Ar 1 (pT 1 ) = 








V 



2sin 2 ^ 
\ Wl 

o -sinl'^ 2 -] 







. / 2pir UJ 2 
sin 

V w i 

2 fpTVU>2 



\ 



Wl ) \ LOl J ) 

satisfies the assumptions of statement (ii) of Theorem [6] because the determinant 



2 sin ' 



w i 

/ 2j)7T W2 



2p7TW2 



. 2 f 2wrw 2 

2 sin ' 



in2 '^r) 



So wc can apply Theorem [5] to system (|f4[) . 

Now £ : R 4 -> M 2 is £,(X, Y, Z, W) = {X, Y). We calculate, when 6^0, the function 

g(X ,Y )=g(a)=d-^-J P 1 JW- 1 (t)Gi(t,x(t,z a ))*J , 

and we obtain the function ^(^Oi^o) 
(17) 

sin (wi r) [26 [ Jfi(r) + A + (tf 2 (r) + F 2 ) sgn [ il^t^ ^ 



2VApTi Jo 
+ (a-b + ^A) (K 3 (t) + F 3 + {K a {t) + Fa) sgn(Xx ,r (r))) 



2/- 



dr. 



and the function ^(^Oj^o) 
(18) 

1 /-p t i 



2x/Apr! 7 



cos («! r) f +2fe f A'i(r) + Fi + (A' 2 (r) + A 2 ) sgn f ^^l Ix ,y (r) 



26 



+ (a - 6 + Va) (A 3 (t) + A 3 + (A' 4 (r) + F 4 ) sgn(X Xo ,y (r)))) 



rfr. 



where the functions of Fi for i = 1,2,3,4 are the ones given in ([5]). Note that -a + d + vA > 0, 

then 

Sgl1 ( ~ — 2b X Xo,Yo( T ) J = sgn(X Xo ,Yo( T )) , 



denoting by A = l/(2VApT 1 ), the function (JTTJ) £7i(-Xo,io) becomes 
-A 



o 



sin (wi r) (A 3 (r) + A 3 ) a - b + VA + 2b{K l {r) +F 1 )dr 



-K I 

Jo 



in (wi r) ((A 4 (r) + F 4 ) (a - 6 + Va) + 26(A 2 (r) + F 2 )) sgn(X Xo ,y (r)) dr, 
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K 



and the function (fT8|) G2{Xq,Yq) becomes 

rpTi . 

cos ( Wl r) (K 3 (t) + F 3 )(a-b+VD)+ 2b{K 1 {r) +F 1 )dr 

jpTi 

K I 

Jo 

Then, by Theorem [6] we have that for every simple zero (Xq,Y *) G V of the system of 
non-linear and non-smooth functions 



i ( Wl t) {[Ki{T) +Fi) (a - b + VD) + 26(^ 2 (r) + F 2 )) sgn(X Xo , y (r)) dr. 



(19) £i(x ,y ) = o , g 2 (X ,Y ) = 0, 

we have a periodic solution (X, Y, Z, W)(t, e) of system (|14l) such that 

pr,Y,Z,W)(0,e)-> (X o *,r o *,0,0) as£^0. 

Note that system (flT)]) is equivalent to system 0, because both equations only differs in a 
non-zero multiplicative constant. 

Going back through the change of coordinates (fT3|) we get a periodic solution (x, y, z, w)(t, e) 
of system (|15|l such that 



( x(r,e) ^ 
J/(i",e) 

z(r,e) 



— {XI cos (wi t) + F * sin (wi r)) 

— (F * cos (ui t) - X * sin (wi r)) 

— (Xq cos (wi r) + Y" * sin (wi r)) 



V 



/ 



y" * cos (wi r) — Xq sin (wi r) 

as e — > 0. 

Consequently we obtain a periodic solution (0i,02)(t,£) of system ([2]) such that 

J -a + b + VA (xx ^ ^ ^ + K „ s . n ^ ^ ^ 
26wi 
1 



(0i,O 2 )(t,e) 



\ 



(Xq cos (wi t) + Y * sin (wi r)) 
(Xq cos (wi r) + Y" * sin (wi r)) 



as e — > 0. Hence Theorem [T] is proved. 



n 



Proof of Theorem [3J This proof is completely analogous to the proof of Theorem Q] 



□ 



3. Proofs of corollaries 

To obtain the expression of the functions given in ((4]) and (0 we have to study the changes 
of sign of the functions Xx a ,Y { T ) an d Z z ° ,w °(t) respectively for t E [0,pTi] and t <G [0,pT2\. 
Note that Xx ,Y (t n ) = for 

Wi V V-Xo 

If XqY) < 0, then i„ G [0,pTi] only for ra = 0, 1, • • ■ ,p+l, and if X Y"o > 0, then t„ G [0,pJi] 
only for n = 1, 2, ■ ■ • ,p + 2. We know that for all t G [t n , *n+i] the function Yx ,y ( r ) nas the 
same sign and different sign for any t € [£„_!,£„], thus the integral can be computed using the 



1:5 



partitions {0, t n , pT\ ; n = 0, 1, • • • ,p + 1} and {0, t n , pT\ ; n = 1, 2, • • • ,p + 2} as the limits of 
integration respectively for XqYq < and XqYq > 0. 

The study of changes of the sign of the function Z z °' Wo (t) for t E [0,pT 2 ] and Z W ^ is 
completely analogous. 



Proof of Corollary [H Firstly, we have to check the Crossing Hypothesis (see Remark [2]) to the 
system ((3|). 

Note that we have four different vector fields defined in four different regions (sec Figure EJ) . 

In the region TZi — {x > and z > 0} we have 



/ 



A', 



V 



\ 



y + 1 

-ax + z + e — ^ h e 2 i?i(i, x, y, z, iu) 

or 



(c 



W7 + 1 



^ bx-bz + e - \- e 2 R 2 {t,x,y,z,w), J 



a 



In the region 1Z 2 = {x < and z > 0} we have 



/ 



A', 



7/ — 1 

—ax + z + e'- — 7. h e 2 i?i(i, x, y, z, iu) 



o- 



(C 



w + 1 9 _ 

V Ox-Oz + e he -R2U, x, y,z,w), I 

or 



In the region 7^.3 = {x < and z < 0} we have 



/ 



X 3 = 



y — 1 „ 

—ax + z + £- — \- e Ri(t, x, y, z, u>) 

1(7 — 1 
I ax — 6z + e - \- e 2 R 2 (t,x,y,z,w), J 



Finally, in the region TZi = {x > and z < 0} we have 



/ 



Xa = 



U 



y + 1 „ 

ax + z + e — = h e Ri(t, x, y, z, w;) 



»' 



id — 1 



^ bx—bz + e — - \-e 2 R 2 (t,x,y,z,w), J 



To study the types of the sets M.^ (see Appendix A), we have to compute Lie derivative of 
the functions g\ and 172 with respect to the vector fields Xi for i = 1,2, 3, 4, i.e. 

(^x t )(9j)(x,y,z,w) = (Vs,-,<*i)(a;,y,z,u>). 
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X = Q 



x% 


*i 


rt) 




^3 ^4 




9i\0) 



z = 



FIGURE 4. Four different vector fields. 



Proceeding with these calculations we have 

(Cxi)(gi)(x,y,z,w) = (£x a ){gi)(x,y,z,w) = y, 

(Cx 2 )(g2)(x,y,z,w) = (Cx 3 )(g2)(x,y,z,w) = w, 
(£x 3 )(gi)(x,y,z,w) = (C Xi )(gi)(x,y,z,w) = y, 
(Cx 1 )(g2)(x,y,z,w) = (Cx 4 )(g 2 )(x,y,z,w) = w. 

Hence we can conclude that in the set 

T = {(x,y,0,0)}\J{(0,0,z,w)}, 

the flow is tangent to the discontinuous set, and in any other point the flow cross the set of 
discontinuity. 

Using the coordinates defined in (fT5|. we have that 



T = 



X,Y, — -X,-Y 

Wl 



u 



X,Y,(3^X,PY 



Observe that the periodic orbits given by Lemma[5]filling the planes {(X, Y, 0, 0)} and {(0, 0, X, Y)}, 
except the origin, do not reach the set T ■ Thus, for |e| > sufficiently small, there exists a 
neighborhood of the planes {(X, Y, 0, 0)} and {(0, 0, X, Y)} such that the orbits cross the set of 
discontinuity. In other words, the Crossing Hypothesis is satisfied. 

Studying the changes of the sign of the function Xx ,y (t) for t £ [0,Ti] we conclude that 
the non-smooth functions (01 and ([5]) are given by 



•Fipf ,r ) 



2VA7T 4(o- 
Y H 



A) 



OJl 



2y/Air 



Ul\ 



Y - 



OJlJl 

4(a + fr + 

W14/1+ 






A) 



* 2 o 



if X Q > 0, 



if X Q < 0, 



lo 



2VAtt 4(o- 
Yq 



•^(-^Oj^o) 



A) 



OJi 



/ Y 2 

"V 1 + s 



2VAtt, 



-Vn 



Wi 



4(a + b+VA) 
I yI 



if X > 0, 



if In < 0, 



This system has all solutions inside a periodic orbit of the unperturbed systems passing 
through 



A7T 



(a + 6 + VA) 



A-7T 



(a + b+VA) 



It is easy to check that this solution are simple. So, by Theorem[T]we have one periodic solution 
of the non-smooth perturbed double pendulum. This completes the proof of the corollary. □ 



Proof of Corollary [7} This proof is completely analogous to the proof of Corollary [H 



□ 



Appendix A: Basic concepts on Filippov systems 

We say that a vector field X : D C R" — > R™ is Piecewise Continuous if the domain D can 
be partitioned in a finite collection of connected, open and disjoint sets Di, i = 1, • ■ ■ , k, such 
that, the vector field X\-^ is continuous for i = 1, ■ • • , k. 

We denote by Sx C 9DiU- • -UdDk the set of points where the vector field X is discontinuous. 
By assumptions, the set Sx has measure zero. 

If Ad C 5a- is a manifold of codimcnsion one, then M. can be decomposed as the union of 
the closure of the regions (see Figure [S]): 

Y, c = {xeM: (Xh)(Yh){x) > 0} ; 

S e = {x e M : (Xh)(x) > e (Yh)(x) < 0} ; 

S s = {xeM: {Xh){x) < c (Yh)(x) > 0} . 






FIGURE 5. Crossing region (S c ), escaping region (E e ) and sliding region (E s ). 
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For p G S e U E s we define the Sliding Vector Field as 

(20) Z s (p) = l _ ((Yh)(Xh)(p) - (Xh)(Yh)(pj) . 

Consider the following equation 

(21) x = X(x), 

where X : D c R™ — > R™ is a piecewise continuous vector field. The local solution of the 
equation (|21j) passing through a point p £ M. is given by the Filippov convention: 

(i) for p 6 E c such that (Xh)(p), (Yh)(p) > and taking the origin of time at p, the 
trajectory is defined as <pz(t,p) = fY(t,p) for t <E I p fl {i < 0} and ipz(t,p) = ipx(t,p) 
for t e / p n {f > 0}. For the case (X/i)(p), (Y7i)(p) < the definition is the same 
reversing time; 

(i) for peS e US s such that Z s (p) ^ 0, (£z(i,p) = <pz a (t,p) for te/ p CR. 
Here (pw denotes the flow of a vector field W. 
For more details about discontinuous differential equation see Filippov's book [J. 

Appendix B: Basic results on averaging theory 

We present the basic result from the averaging theory that we shall need for proving the 
main results of this paper. 

We consider the problem of the bifurcation of T-periodic solutions from differential systems 
of the form 

(22) i(t) = G (t, x) + eG^t, x) + e 2 G 2 (t, x, e), 

with e = to e ^ sufficiently small. Here the functions Gq,Gi :1xS1-> M™ and G2 '■ 
Kx O x (— £o,£o) — > K™ are C 2 functions, T-periodic in the first variable, and f2 is an open 
subset of R n . The main assumption is that the unperturbed system 

(23) x(t)=Go(t,x), 

has a submanifold of periodic solutions. A solution of this problem is given using the averaging 
theory. 

Let x(t, z, e) be the solution of the system ([23]) such that x(0,z,e) = z. We write the 
linearization of the unperturbed system along a periodic solution x(i, z, 0) as 

(24) y = D x G o (t,x(M,0))y. 

In what follows we denote by M z (t) some fundamental matrix of the linear differential sys- 
tem (f24|) . and by £ : K fc x M" _fe -> R k the projection of R" onto its first fc coordinates; i.e. 
£(cci,...,a; n ) = (xi,...,x fc ). 

We assume that there exists a fc-dimensional submanifold Z oi ft filled with T-periodic 
solutions of (|23[) . Then an answer to the problem of bifurcation of T-periodic solutions from 
the periodic solutions contained in Z for system (f22j) is given in the following result. 

Theorem 6. Let V be an open and bounded subset of~R k , and let (3 : C\(V) — ¥ M"~ fc be a C 2 
function. We assume that 

(i) Z = {z a = (a, /3(a)) , a € C1(V)} C fl and that for each z Q € Z the solution x(i,z Q ) 

of ((23)) is T-periodic; 
(ii) /or eac/i z Q G -Z f/iere is a fundamental matrix M Za (t) of (|24[) suc/i t/iai f/ie matrix 
M~ 1 (0) — M^iT) has in the upper right corner the k x (n — k) zero matrix, and in 
the lower right corner a (n — k) x (n — k) matrix A Q with det(A a ) 7^ 0. 
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We consider the function Q : C1(V) 

(25) Q{a) = t(^J M-^tJdCt.x^aa))* 

If there exists a <E V with G(a) — and det ((dQ/da) (a)) 7^0, i/ien i/iere is a T -periodic 
solution <p(t,e) of system (|22|) s«c/i i/iai yj(0, e) — > z a as e — > 0. 

Theorem |6] goes back to Malkin [6] and Roseau [7] , for a shorter proof see [3] . 
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